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By Harvey G. McComb, Jr. 
SUMMARY 


Equations are derived for the stress distributions caused by three 
types of loading on infinitely long, circular, semimonocoque cylinders 
with flexible rings. These external loads can be said to lie in the 
sheet-stringer surface of the cylinder; that is, they have no radial 
components. The results, given as formulas for the stringer loads and 
shear flows in the shell, are of use in the stress analysis of circular 
semimonocoque cylinders with cutouts. In order to facilitate this appli- 
cation, the formulas can be used to construct tables of influence coef- 
ficients giving stringer loads and. shear flows in the neighborhood of 
each particular external load due to & unit magnitude of that load. 


INTRODUCTION 


An airplane of semimonocoque construction usually has various 
openings for doors, cargo hatches, landing-gear mechanisms, and other 
purposes. The structural design of the airplane near such openings 
requires knowledge of the stress distribution about cutouts in semi- 
monocoque shells. Stress analysis of semimonocoque structures with 
cutouts is beyond the scope of engineering beam theory because stress 
concentrations occur in the neighborhood of the cutouts. it is for 
this reason that more powerful methods of analysis are sought. 


A method for the stress analysis of sheet-stringer panels with cut- 
outs has been presented in reference 1. The method is based on the 
idea introduced by Cicala (ref. 2) that the effect of the cutout can be 
reproduced by superposing certain perturbation stress states on the 
stresses which would occur in the structure without a cutout. Three 
types of “unit perturbation solutions" were obtained, and it was shown 
how these solutions could be used to handle cutouts in pee eee 
panels under axial loads. 
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It is evident that this approach can be extended to the analysis of 
stresses around cutouts in circular semimonocoque cylinders. In order = 
to make the extension, however, perturbation solutions for this type of 
structure must be developed. Accordingly, the purpose-of this report is 
to derive the perturbation solutions for an infinitely long circular 
cylinder of semimonocoque construction. These solutions are the stress 
distributions due to three types of loads analogous to those discussed 
in reference 1. The particular applied perturbation loads for this 
structure are shown in figure l. In the analysis for the stresses due 
to these loads, bending flexibility of the rings in their own planes is 
taken into consideration. 


The use of these perturbation solutions in the stress analysis of 
circular semimonocoque cylinders with cutouts is described in detail in 
reference 3. 


SYMBOLS 
A cross-sectional area of stringers 
= SES - 1 + cos n5 
fati 
ge Avifn(i) (n > 2) : 
2L sin ae 
2 
bt e 
By = 3B6° + 2(1 - cos nb) 
b arc distance between stringers, Rẹ 
f (i 
bai = - DA (n > 2) 
2L sin T 


ARÓ 


C = = 


IL?b 


Ep ei 


Ze 
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| 2(Bn - 1) 
yc 

TE E 

(zm + n)? Kem + n)? - | 
dn coefficient in trigonometric series for Do: 
E Young's modulus of elasticity 
F; j tangential force on ring i uniformly distributed between 

stringer j and stringer j + L 

EE coefficient in trigonometric series for stringer loads 
G shear modulus of elasticity 


H,(n,%) = = Dan cos(rm + n)ó 


T==0 
Ho(n,g) = > (-1) Dan sin(rm + n)$ 
I-A) 
I moment of inertia of ring cross section 


i longitudinal index, indicating stations (where rings are 
located) and bays 


Jy has the value l when r is an integer and has the value 0 
when r is not an integer 


j circumferential index, indicating stringers and panel rows 
k,1,r,8 integers 

L distance between rings 

M1,9) ` bending moment in ring i 

m total number of stringers in the shell, m 23 


n index of the terms in a trigonometric series 
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P external concentrated force in the longitudinal direction 
applied to a stringer at its intersection with a ring 


Pij stringer load in stringer j at ring station i 

Q external shearing force per unit EES SECH about a 
shear panel | T 

dij shear flow in shear panel (i,j) 

R radius to middle surface of sheet i 

S external force in the longitudinal direction uniformly 


distributed along that portion of a stringer which lies 
between adjacent rings 


y-—-oo 
T(1,9) thrust in ring i 
t thickness of sheet 
U total stress energy 
V(i,@) transverse shear in ring i 
O1n»92pn 
arbitrary constants 
O37 An 
2 
LA 3 BS e 
sn È 
Bn = 3 + : 
l2C84 E 
ba B5. 
S sin? 2» 
Yn 7 -2 + | 
leCcSa 
Aii second central difference in the i-dírection or longitudinal 


direction, that is, Asye(i) = gli + 1) - 2g(i) + g(i - 1) 


i 


ku, 


Cre 


d: 


a, 
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5 central angle between stringers, 2/m 
Org Kronecker delta, takes the value 1 when r = 8 and takes 
the value O when r #5 
t, = te "n 
Andon quantities defined immediately following equation (21) 
g angular coordinate for rings 
Së p l Bn + 1 S 2 
= -1|Bn - 
Xn = 5 cosl E - (nt) - ?n (Dy > 1) 
= l eosnl Pa cl (ea + 2) e (D, < 1) 
ud —$ "I$, 7 7n Dn 
- ed / + 1\¢ 
Yn = 4 cosh DE + ("zt ) "a 


BASIC ASSUMPTIONS 


A typical portion of a structure of the type discussed in this 
report is illustrated in figure 2. It consists of a thin-walled circu- 
ler cylinder stiffened by stringers in the longitudinal direction and 
by rings in the circumferential direction. The structure is infinite 
in length. 


The analysis in this report is based on the following assumptions 
regarding. the properties of the structure: 


(a) The stringers are uniform and equally spaced around the shell, 
and the sheet is of constant thickness. 


(b) The stringers carry only direct stress. The sheet takes only 
shear stress which is constant within each bay; thus, the stringer 
stresses vary linearly between rings. 


(c) The rings have a finite bending stiffness in their own planes, 
but they do not restrain longitudinal displacements of the stringers. 
The ring bending is inextensional. 
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(a) The difference between the radius of the neutral axis of the 
rings and the radius of the middle surface of the sheet is negligible. co 


ANALYTICAL APPROACH 


Equations for the stress distributions arising from the three per- 
turbation loads illustrated in figure 1 are derived in this report. The 
perturbation load shown in figure l(a} is a concentrated force in the 
longitudinal direction applied to one stringer of- the shell at-a ring 
station. In figure 1(b) the perturbation load shown is a force in the 
longitudinal direction uniformly distributed along the portion of one 
stringer of the shell between adjacent rings. Finally, the perturbe- 
tion load shown in figure l(c) consists of forces uniformly distributed 
on the portions of stringers and rings which border one shear panel of 
the shell, the forces acting in such a manner as to cause pure shear 
in that panel. These three loads will be called, respectively, the 
concentrated, the distributed, and the shear. perturbation load. 


A detailed explanation of the use of these perturbation loads in 
the stress analysis of circular semimonocoque cylinders with cutouts is 
presented in reference 3. Very briefly, the idea is as follows: The 
cylinder is analyzed by any applicable method with the assumption that 
the cutout does not exist. Then a system of perturbation loads is 
applied to the shell without the cutout, and these loads are adjusted » 
in magnitude to force the boundary conditions at the cutout to be | 
satisfied. Stresses are obtained by superposing the perturbation stress 
distributions on the stresses which would exist in the structure without. D 
a cutout. The method is sufficiently versatile to handle the effects of 
reinforcement about the cutout. 


The perturbation solutions will be obtained by use of the principle 
of minimum complementary energy. This principle states that, among all 
possible stress distributions in the structure which satisfy equilibrium 
and the boundary conditions on stress, the distribution that most nearly 
satisfies compatibility is the one which minimizes the complementary 
energy a" where 


= Internal energy - Work done by surface stresses acting through the 
prescribed displacements 


(1) 


Since displacements are not prescribed anywhere on the structure, the 
second term on the right-hand side of equation (1) is omitted. The com- 
plementary energy becomes the internal ey) or stress energy of the 
structure. | = EUR 


Wi 


1 Gë 
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In writing the equation for the stress energy, the following factors 
are considered: the energy of axial distortion of the stringers, the 
shear energy in the sheet, and the bending energy of the rings in their 
own planes. In figure 2 are shown the index systems used to designate 
stringers, rings, bays, and panel rows. Note that the intersection of 
ring i and stringer j occurs at the lower left-hand corner of shear 
panel (i,j). Each of the perturbation loads is shown in its positive 
sense in figure 1. Stringer loads are taken as positive in tension. 
Shear flows are positive as shown in figure 2. Ring bending moments, 
shear, and thrusts are placed on the ring element in figure 3 in the 
positive sense. With this index notation, the stress energy in the 
structure can be expressed as 


oo 


m= | | | 
z d 2 2 RSL 2 
U a" Za, + PijPisi,j + P141,32) * 2t Td + 


dE 

>f zz; M(i,p)ag (2) 
Zeep 2æÆI 

Where the integration over the length of a stringer between adjacent rings 

has been carried out. 


In the analysis to follow, stringer loads are expressed in the form 
of a finite trigonometric series. Then, by using the equations of statics, 
the shear flows and ring bending moments are written in terms of the 
coefficients of this trigonometric series. The expression for stress 
energy, equation (2), is minimized with respect to the coefficients of 
the trigonometric series for stringer loads; then, the expressions for 
the stringer loads, shear flows, and ring bending moments are substituted 
into the resulting equation. This process yields a fourth-order finite- 
difference equation which can be solved for these trigonometric coeffi- 
cients. The solution is then substituted back into the original expres- 
Sions for stringer loads, shear flows, and ring moments to yield the 
desired distributions. 


For convenience in application, the significant equations are col- 
lected in appendix A. 
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PERTURBATION LOAD SOLUTIONS 


Concentrated Perturbation Load 


ression for_stringer loads.- The concentrated perturbation load 
is shown in figure l(a); let P represent the magnitude of this load. 
Since the structure is uniform and infinitely long, half of the load 
goes into the portion of the structure to the right-of the ring station 
where the load is applied (ring i = 0), and half goes to the left of 
this ring station. Therefore, it can be seen from figure l(a) that, 
because of symmetry, 


Pij = “P-i,3 | (121) (3a) 
dig ed coq d (120) (3b) 
M(1,8) = -M(-1,9) (120) (3c) 


Consider the right half of the structure, including the ring at 
i = 0. The concentrated perturbation load gives rise to stringer loads 


which are circumferentially symmetric about-stringer j = O (see fig. 1(a)). 


Thus the stringer load distribution can be represented by a series of the 
form i 


m m=]. 
2 7$ 
Pij = > £,,(1)cos njo l (4) 
n=0 | | 
Moy MeL 
e 2 
where the notation > means that the summation is carried over n 
=O 
from n-0 to n= 5 1f m is even and to ne 22 if m is odd. 


Evaluation of: fo(i), f,(1), and fy(0).- Suppose that equation (+) 
is multiplied by cos 2j& and summed over j from O to m- 1. This 
procedure yields l | ES 


p13 cos ij5 = SN fn(i) M cos nj8 cos 1j5 
J= n=0 j=0 


rt 
I 


di 
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The sum over j on the right-hand side is, for OS pn 


lA 
MIE 
d 
O 
WA 
C4 
WA 
hole 


H 
© 
om 
CJ 
“tH. 
3 
c 


E cos DAD cos Łjò 
j=0 


In z= 
E + no t 55,2) (1 = n) 


Thus the coefficients of the trigonometric series in equation (UI are 


m= 
2 $ 
fali) = ————————— p,4 COS njb (5) 
el + no + 5,2) = 


It is desirable first of all to evaluate those f,(i)'s which can 
be found from consideration of the boundary conditions and of the over- 
all equilibrium of the cylinder. Consider the equations of statics for 
the cylinder as a whole. Satisfaction of equilibrium in the longitudinal 
direction requires that the sum of the stringer loads at any ring sta- 
tion i mist equal one-half of the applied load P. This condition is 
expressed as 


For n - 0, equation (5) gives 


fo(1) = i 2 Pig = oe (6) 


Moment equilibrium gives two equations, one of which is automatically 
satisfied because of the symmetry of the stringer load distribution around 
the cylinder. The other moment equation is 


-1 
S Dua R cos Jo = PR 
dei 
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For n = 1, equation (5) is i 


£1(1) =£ a p,4 cos JO T | (7) 


On substituting the values of fo(i) and fı(i) given in equations (6) 
and (7), respectively, into equation (4), there results 
M or IER) 
P P 


: m 
Dij = ga t g 008 jò + > f a(i)cos nj (8) 


Consider now the boundary condition at ring station i = 0. The 
stringer loads here are 


"d 


Poj = 5 Do 


and substitution of this expression into equation (5) yields 


fy(0) = ucibus (035232) (9) 


ns + 


Ihe equations of equilibrium and the boundary condition at i = O 
have been used to obtain certain of the coefficients of the trigonometric 
series for stringer loads. The remainder of the coefficients are found 
by use of the principle of minimum complementary energy, and this is the 
next step in the solution. 


ressions for shear flows and ri bendt moments.- In order to 
use the principle of minimum complementary energy, the shear flows and 
ring bending moments mist be found in terms of the trigonometric coeffi- 
cients fp(i). Shear flows are determined by the consideration of-the 
equations of statics of a portion of any stringer j between two adjacent 
rings i and i+ 1. The forces on this free body are shown in the sketch 


below: 
Oy at 


— Pun 
Pig —— Pi, 
Wë, EE 
dr jet 


L 
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Equilibrium of these forces requires that 
Pi+l,j 7 Pij + (a3 - 44, 3-1)b -0 (10) 


Substitution of equation (8) into equation (10) yields 


TT 
1 
34-8,41-7-i pa [f (11) - f4(1)]eos afp (11) 


In order to find Q4 j? replace j with a dummy index k and sum both 
sides of this equation over k from kK=1 to k= j; that is, write 


op Met 


I 
) 2 j | 
2 (ax. J di ke) == = gës [en (1+1) - SCH ex cos nk5 


MIB 


When the indicated summations over k have been carried out, the fol- 
lowing equation is obtained: 


1 
—— ee N Kä) - sti SE «i 


L 


The term on can be found from the condition that the total torque on 
the section is zero. The resulting expression for shear flows is 


m m-1l 

2 OF "5 fp(itl) - fali) , ( De Se 
daa = - ——— gin n += 
ij n-2 2L sin Si e 


Bending moments are caused in each ring by a tangential loading 
which develops because of the difference in shear flow in the sheet on 
either side of the ring. The tangential load on ring i has the value 


m oy Mel 


2 2 e (li 
SEES es ey er 2 Alifnit/ nt sin sl + 2 (15) 
^ n=2 eL sin = e 
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In appendix B, this load is applied to a circular ring and the bending 
moment in the ring is derived. This procedure results in the following 
moment in ring i (see eq. (B9)): 


M(i,g) = - es Rm Aifa(1)(m, 9) — (11.) 
n- 


where 


(rm + n)ó 
> cos 
r--9 (rm + n)? Ven A n)? g 1] 


The sign convention for the moment is illustrated in figure 5; the con- 
vention for measuring angle Ø is shown in figure h(a). 


Energy analysis.- Now the stringer loads, shear flows, and r 
bending moments have been expressed in terms of the coefficients fa i). 
The stringer loads are given in equation (8), the shear flows in equa- 
tion (12), and the bending moments in equation (1+). These equations 
are used in the minimization of the stress energy of the cylinder with 
respect to f,(1). 


By virtue of the symmetry properties of this problem expressed. in 


equations (5a), (3b), and (3c), the energy in the structure to the left - 


of ring i = O is the same as the energy to the Eer of ring i = 0. 
Thus, Eoun (2) becomes 


al, S on. RSL 
u=2>_ (es t PijPi+l,j * 2141,32) ` Se vil j 


i=0 


dëi SIS M (1,9)ad 


Note that M(0,@) is identically zero because there is no difference in 
shear flow across ring i = O and, therefore, no tangential load acts on 
this ring. 


Vi 


af 
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Minimization of the stress energy with respect to fp(i) results 
in the following equation: 


aU 


Ofy (1) CR 


9p, j Së SCT 
1 
Ee ol GER + 1p: 3 t P1-1,3)s T nti) ate Ap Së KS 


SN oM(i+1,¢) 3M(1,9) , 
d1-1,j Se 2) mos. L Es Hag) Ofn(i) * MSS T n(i) 


OM(i-1,9) 
Ofal i) 


The coefficients fg(i) and fı(i) are known already for all values of i, 
and f,(0) is known for O<n< 2. Thus equation (15) needs only to be 


M(i-1,¢) ag (15) 


considered for i Z 1 and n22. The expressions for the stringer loads , 


shear flows, and ring bending moments are substituted into equation (15). 
Then the following definite sums are needed (these can be obtained by 
the procedure outlined in ref. 4): 

m=] 


j-0 


cos DAP iEn (o <n<m) (16) 


and for the integers n and 2 restricted to the range Lens z and. 
1$7$- 


= p? 
m= 1 
> cos 1j5 cos njB = 0 (i d n) 
Ze? 
(17) 
E a(t + dÉ (2 =n) 
and 
2 ein 1 + Bp sin aG ec (2 #n) - 
| i 


jl 
VIE 
Z^ X 
T 
o 
p 
D 
bU d 
gen 
| 
B 
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The following definite integral, which is derived in appendix C, is also 


needed: ^. 
2x 
T H1(1,2)831(n,9)ad = o (1 f n) 
O 
(19) 
= el 4-5 m (2 = n) 
Has 
where 
CO oo 
2 1 
Sn = > Din = > mo hh, es TB 
EE r=- (rm + n) rm + n)* - " 
and where n and 2 are restricted to 2$n$7 and 2s LSS A 


closed form of S, is presented in appendix C but the series form con- 
verges so rapidly that it is usually more convenient than the closed 
form for use in calculations. 


After substitution of the expressions for stringer loads, shear 
flows, and ring moments into equation (15), the use of these definite 
sums (16), (17), and (18), and definite integral (19) results in the Š 
following equations which express the conđition of minimum stress 
energy: 


For i= 1, 
fn(3) + 2ypfn(2) + (2Bn - 1)fn(1) + 2(yn + 1)£,(0) = 0 (20a) 


and, for i 2 2, © 


fQ(i-2) + Q7pfy (itl) + 984f4(1) + 2y4f4(1-1) + f£,(1-2) = 0 (20b) 


where 


by 
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Solution of finite-difference equation.- Equation (20b) is & fourth- 
order finite-difference equation with constant coefficients. (Note that 


the symbol i represents the index of the rings and bays and should not 
be confused with the usual notation for  /-1 which never appears in 
this report.) Equation (20b) corresponds exactly with equation (2h) of 
reference 5. The general solution is presented on pages 25 to 26 of 
reference 5 and on pages 20 and 29 of reference 6. It may be written as 


fn(i) = (ae a) ERO T conden (1) | + 
(te¥n)* [oann + on (1) | (n 2 2) (21) 


where the upper sign is used when 7, « O0 and the lower sign when 
yn ^ O0. The A's are as follows: 


2 es Í 
For Ea E EN 
Yn* 


Ajg(i) = cos ix, 


Aen (i) = sin iX, 
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- m 


where 
X4 = Ł cos"l Pn- t pa ta 2 
n=3 78 “fC a) 7% 
for D,< 1, 
Ain(i) = cosh ix, 
Ap (1) = sinh 1X4 
where 


e | 
— 1l -l|Pn 7 1 Pn + J 2 
Xn = z cosh mco = Yn 


and, for Dy, = 1, 


ll 
H 


Ain 


ai 


Aon 


ln the inverse trigonometric and hyperbolic functions, the principal 
values are used. The argument yy of the exponential function is given 
by the positive branch of. 


2 
D ai í +1 | 


At a large longitudinal distance-from the applied load, the stringer 
loads should approach the elementary distribution given by the first two 
terms of equation (8); consequently, for n 22, fg(1)—90 as 1i—»o. 
The first term on the right-hand side of equation (21) satisfies this 
condition; however, the second term does not and, hence, must be omitted. 
Ihe solutions, then, that are compatible with the boundary conditions at 
infinity are: ` | ~~ | | ] 


1 


Cafi) = [nana + aanholt) ] (n 22) (22) 


Hj 
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where 
bn = te" Yn 


Now the arbitrary constants aj, and ado, are determined. The first 
aj, is obtained immediately. Substitution of i = O into equation (22) 
and use of equation (9) to evaluate f,(0) yields 


fQ(0) = ain (n22) (23) 


ee: ENEN 
m(1 + 55,2) 


Substitution of equations (23) and (22) into the boundary equation (20a) 
yields l 


- 614 + 207, + 1) p 
l+ 
85, n( "oP 


where 


Osn = E OG) + Sab "Aen (2) + (284 - 1) nAsn(1) (s = 1,2) 


The solution for the concentrated perturbation load is now complete 
since the coefficients f,(i) are completely defined and may be substituted 
into equation (8) to give the stringer loads. The shear flows can be 
found from equation (12); however, once the stringer loads are known, 
it is simpler to calculate the shear flows by the use of the equations 
of statics. Because of symmetry, the shear flows in shear panels adjacent 
to stringer A = O are given by 


PLO 7" P14,1,Q 
= mA ——ÓMÀMÁ 
WO 7 781,-1 Es 


All the other shear flows can be found by the use of equation (10). If 
desired, the moment distribution in the rings can be computed from 
equation (14) and the thrust and transverse shear in the rings can be 
found from the formulas given in appendix B. 
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Distributed Perturbation Load 


ression for stringer loads.- The distributed perturbation load 
is shown in figure 1(5); let S represent the magnitude of the total 
force distributed along stringer A = O between rings i = O and i = 1l. 
From figure l(b) it is seen tbat 


Pig = 7P-i141,4 l (i 2 1) (24a) 
M(i,@) = -M(-i+1,øØ) (4 2 1) (cl 


At ring station i = 1 and to the right of this ring, the stringer 
loads can be represented by a trigonometric series of exactly the same 
form as equation (kh) 


m M-L 
e 


2 
Pij = 2. fy(ijeos nJ5 (25) 


n=O 


except now i Z 1, and the coefficients fy(i) are different from those 
obtained for the preceding case of the concentrated load. 


Evaluation of fo(i) and £f4(i).- As in the preceding case, the 


first two coefficients fo(i) and fi1(i1) can be obtained from the equa- 


tions of statics, and the results are the same as before. Equation (25) 
becomes to " d . - E See un 


m m-1l 
n 5" 
S 
Pig = oe + P cos jò + gd fy(ijeos nj (121) (26) 


With the concentrated perturbation load, all the coefficients f,(0) 
were easily found because the stringer load distribution at ring sta- 
tion i = O was known. Here no such distribution is known. In order to 
determine the boundary condition at bay i = 0, the effect of the dis- 
tributed perturbation load on the equilibrium of portions of stringers 
in this bay must be investigated. 


Expressions for shear flows and ring bending moments.- Away from 


bay í = O the shear flows and ring bending moments are of the same form 


(d 


«t 
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as for the concentrated load. The following expression for the shear 
flows is obtained by use of equation (10): 


M or m 
2 SCH fa(i-l) - fQ(i) 


n-e 2L sin Gi 


WV 


q13 7 sin d A Ss (4-21) (27) 


The ring bending moments are obtained in appendix B as 


m m-1l 
= or == 
D r 


2 
M(i,g) = - Bn ROD Ay fn(1) (np) (a 
n- 


iv 


2) (28) 


where 


CO 


Hy(n,g) = A ee e 
SE (rm + 22 | (rm +n)*- 1| 


The applied force in bay i = O may be written as S504* Consider, 
now, the equilibrium of a portion of any stringer j between ring i = 0 
and ring i = 1l. The forces on this free body are shown in the sketch 
below: 


Po; M—— RB > P4; 


Eel 


Equilibrium of these forces requires that 


P313 7 Poj + (203 7 ao, 3-1)P - oj 70 
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Because of the antisymmetry property expressed in equation (218), the 
equilibrium equation becomes 


2713 + (303 - We? EE (29) 


It is convenient, now, to expand the Kronecker delta Do: in a 
finite trigonometric series, | ` 


503 = > dy cos njb (30) 
n=O 


Multiplying through by cos 1j6 and summing over j from O to m- l, 
yields the trigonometric coefficients dy» The result is 


dy LÁ A (31) 


n» 


Substitution of the expression for stringer loads, equation (26), and 
the trigonometric expansion for 503; equation (30), into the equilibrium 
equation (29) yields 


037 4,4317 p 2. (B - 2fn(1)] cos ngo 


n-2 


In order to find doj? this equation can be treated in the same manner 
as equation (11); that is, replace 3 by a dunmy index k, sum from 
k=1 to k= j, and then use the condition that-the total torque on 
& cross section in bay i = O must be zero. This procedure. results in 
the following expression for the shear flows in bay i = 0: 


zl 

Bor e i Sd, - fy(1) 2 

doy = ds Ne x sin dÉ + Ls (32) 
= n E 


Ihe expression for tbe bending moment in rings i = l and i = O is 
yet to be found, as this expression differs from that for the moment- in 
the rest of the rings given in equation (28). The moment in ring i = 0 
is the same in magnitude as that in ring 1 = 1 but opposite in sign. 


^n 
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The tangential loading on ring i = 1 is given by 


m or m-i 

2 FS 2 f,(2) - 3£,(1) + Sdy : ( S iy 
q ve = = ———— GIN EL el - 
lj "o SE 2L sin BS e 


2 


By analogy with equations (13) and (14), then, the bending moment in 
ring i = 1 can be written as 


m=] 
a. 


2 
M(1,9) = - 2. Zb - 3fn(1) + Sàg [my (m, ø) (33) 


All the stringer loads, shear flows, and ring bending moments have 
now been expressed in terms of the coefficients fu(i). The stringer 
loads are given in equation (26), the shear flows in equations (27) and 
(52), and the ring moments in equations (28) and (433). It remains to 
substitute these expressions into the equation obtained from minimization 
of the stress energy of the cylinder with respect to f,(1i). 


Energy analysis.- By virtue of the symmetry properties in this 
problem, given in equations (24), the energy in the structure to the 
right of bay i = O equals the energy to the left of this bay.  Equa- 
tion (2) for the stress energy cen be written 


m= L ce M- 
Me 2 (dis Ps? + Gi doj d GER Si + P1jP141,j * Pil j 5" 


i=l 


Eg E (2D Là Witz, dia (34) 


act " 4=1 


Minimization of the stress energy with respect to fnf{i) results in the 
following equations: 
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3U | | E » 
Of (1) diu 
- s 3 2c òq] j Ze 
j ai pij * sj) LE) t Sat Se CL) * doj EE 
ex aM(2 9), oM(1,¢ 
R , 
J = kee, M(1,8)— 1 er d (35) 
and 
OU 
9fg(1) 


+ BBL 905 
Pu w E -ü p SÉ 8r.) 


SML, d) l AM(i, Ø) 
d1-1,j Sr una). I ET ca e Ses) T Wil + 


deeg ag (1 2 2) (36) 


Note that equation (56) is the same as equation (15), except here it is 
valid only for i 2 2. 


The stringer loads, shear flows, and ring moments are substituted 
into equations (35) and (36), and then the definite sums and definite 
integral derived in the preceding section are used to simplify these 
equations. After simplification, the following equations result: 


For i= 1l, 


- hyn - 2 
fn(3) + (27n - 1)fp(2) + 2(Bn - 7n)fn(1) = sa E (37a) 
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for i = 2, 


fn(h) + 2y4£4(5) + Bpfnl2) + (27, - 1)fn(1) = -Sdn (37b ) 


and, for i 2 3, 


Py(it2) + 2ygfg (i41) + 2Bnfn(i) + 27pfn(i-1) + fp(i-2) = 0 (37c) 
Solution of finite-difference equation.- Equation (37c) is the same 
as equation (20b); therefore, the solution to equation (37c) is 
fnli) = bn [einAin(4) GI (n22) (58) 
n n |C1n^in + Gorton nce 5 


which is the same as equation (22) except for the values of the arbitrary 
constants ay, and app- These constants are found by the substitution of 
the solution (38) into equations (37a) and (51b). This procedure yields 
two simultaneous algebraic equations in ay and apy, and their solu- 
tion gives 


Bn - byn - 2 zu 
2n S 2n - 
Lin = 
Anlon - Tinton WE m) 
n5 
Bn - Hyn - 2 
au tig ———— 28 
apn Ss ————— 
Dia ~ Finjen si + 55,2) 
; "S 


where dg, the coefficient in the trigonometric series for the Kronecker 
delta 993» has been replaced by its value as given in equation (31), 
and where the Q's and I's are given by 


Men = to^ Aen (3) * (2x4 - UE Aen(2) + 2(Bn - 7n)tn^en(1) (s = 1,2) 
(592) 
len = Eo hg) t y bn Aent 2) + 28ntn“Agn(2) + 


(27n - 1)6 Asn(2) (s=1,2) (39) 
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The coefficients fp(i) are now defined for the distributed perturbation 
load and may be substituted into equation (26) to give the stringer loads. 
The shear flows can be found from equations (27) and (32), but, again, 
once the stringer loads are known shear flows can easily be found by 

use of the equations of statics. The shear flow in the panels adjacent 
to stringer J = O can be found by considering symmetry: 


In bay i = 0 


S - 2p 
doo E "Hy sec ox 


and, outside of bay i = O, 


Pio 7 Pia1,0 
2L 


ug 94. (1 2 1) 


The other shear flows are found from equation (10), as before. If 
desired, the ring moments can be obtained from equations (28) and (33) 
and the ring thrust and transverse shear can be-found from the equa- 
tions given in appendix B. 


Shear Perturbation Load 


ression for stringer loads.- The shear perturbation load is shown 
in figure Lei, The magnitude of the load per unit length applied along 
the stringers and rings bordering shear panel (0,0) will be represented 
by Q. From figure l(c) itis seen that the longitudinal symmetry prop- 
erties in this case are the same as those for the case of the distributed 


perturbation load given by equations (2L). 


The shear perturbation load is self-equilibrating and gives rise 
to stringer loads which are antisymmetric about panel row j = O. For 
i 21, the stringer loads may be represented by 


m mel 
= or == 
2 


2 
Daa = f,(1)sin n(j - Łe (40) 
SET wech? 


Where the coefficients f,(1) are different from those in the two pre- 
ceding cases. The term corresponding to n= 1 vanishes because it- 
represents an elementary bending stringer-load distribution, and the 
shear perturbation load does not require this distribution for overall 
equilibrium. 


k; 


r1 
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ressions for shear flows and ring bending moments.- None of the 
coefficients Pa (i) in the trigonometric series Zico) can be found from 
the equations of statics. Furthermore, the boundary condition at | 
bay i 2.0 must be determined from a consideration of the effect that the 
shear perturbation load has on the equilibrium of the portions of stringers 
in bay i = O and on the bending moment in the rings bounding this bay. 
Thus the energy approach must be used immediately and the first step in 
this approach is to write the shear flows and ring moments in terms of 
f,(i), the coefficients of the trigonometric series for the stringer 
loads, equation (40). 


Outside of bay 1 = 0, the satisfaction of the equations of statics 
for the portions of stringers between adjacent rings yields equation (10), 
the same as in the two preceding cases. Substituting equation (40) for 
the stringer loads into the equilibrium equation (10) and following the 
Same procedure used to obtain equation (12) yields the expression for 
the shear flows due to the shear perturbation load: 


In m-l 
2— 2 f (1+1) - fali) 
qij = ————————— cos nj5 (i 21) (41) 
n= 2L sin = 


The tangential loadings on the rings to the right of ring i = 1 are 


m m= 1 


2 SS 2 Aiifn(i) 


n=2 eL sin © 


In appendix B this load is applied to a circular ring and the following 
expression for the moment in the ring is obtained (see eq. (B13)): 


cos nj 


4:3 7 &1-1,3 7 


o or SC l 
Mp) = TS" BB aur. RS) (222) (e) 
n=2 E 


where 


co 3 n)g 
(n,9) = > -1)f .ooLB8in(rm + n) ` 
a poss (rm + n)? Krm t n)? - i| 


The convention for measuring angle Ø here is a little different than 
before and is illustrated in figure }(b). 
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Now, the shear flows in bay i = O and the bending moments in the 
rings bordering bay i = O must be found. Consider the shear flows in 
this central bay. The shear perturbation loading applied at bay i = O 
may be written Poz. Then the forces on the portion of any stringer 3 


between ring i = O and ring i = 1 are as shown in the sketch below: 


L 
P ds 
Goal, "4 
zb Bn Di: 
Es 
L N 
do, 3-1 v 
Bo 3-1L 
L 


Equilibrium of these forces requires that = 
Ei: e Poj + (26; es do, 4-1) T ao, 4.1 a Soa =O 


Because of the antisymmetry property, equation (24a), the equation of 
equilibrium becomes 2 


s F Loo: S 30, 3-1)L S q(5o, 3-1 - doj) = O (43) 


The substitution of the stringer loads, equation (40), into the equilibrium 
equation (43), and the introduction of the trigonometric expansion for the 
Kronecker delta Bou: equation (30), yields the following equation: 


m m- 1. 
2 or -57 ; 
A 1 
193 7 90,j-1 = “T ) 2fp(l)sin sl, - 2: - Qài cos(j - 1)8 - cos J E 
n-2 
m nn m-l 


: >. Qi, [cos n(J - 1)8 - cos nj] 
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Now doj can be found by replacing j with a dummy index k, summing 


over k from k=1 to k= j, and using the condition that the torque 
on & cross section within bay i = O balances the applied torque. This 
procedure results in the following equation for the shear flow in the 
central bay: 


-2 € fall) 
doj = Wo + Mı cos jè + PN + sde nj5 (L) 


Consider the bending moment in rings i = l and i = O. The moment 
in ring i = O is identical in magnitude to the moment in ring i = 1 but 
of opposite Sign. The tangential loading per unit are length on dens 
portion of ring i = 1 between stringer j and stringer j + lis illus- 
trated in the sketch: 


Bo 
Ve 
T — `~ = =- § stringer J+ 1 
D Go: dj 


bk stringer j 


Adding these tangential loads and introducing the series expansions for 
304 and 215 and for Go: yields the total load per unit arc length 


on ring i= 1: 


2 "2 ciel. 3£n(1) 

fn(2) - 3fn(1 
Gis e + Co4 = — cos nj5 
1j 7 703 " ER 2L sin È 


2 


By analogy with equations (13) and (14) the bending moment in ring i = 1 
is 


5 


Wg) = SGT Rot rato) - 5£2(2)]ge(n,9) pa 
n= 
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Expressions for stringer loads, shear flows, and ring moments have 
been written in terms of the coefficients fp(i). The stringer loads are 
given in equation (40), the shear flows in equations (41) and (44), anda 
ring moments in equations (42) and (45). These expressions are ready 
to be substituted into the equation which results from minimizing the 
stress energy with respect to fali). 


Energy analysis.- Because the longitudinal symmetry relations which 
exist for the distributed perturbation load, equations (24), also exist in 
the case of the shear perturbation load, the stress-energy expression used 
in the distributed-load problem can be used here. The expressions obtained 
on minimizing this stress energy, equations (35) and (36), are also appli- 
cable here. So the stringer loads, shear flows, and ring moments just 
derived are substituted into equations (55) and (36). At this stage in 
the two preceding cases, certain definite sums and a definite integral 
were introduced to simplify the equations. A similar procedure is 
followed here. - M | | 


The definite sums which are of interest are. 
mM- JL 
> sin sl, - Le = 0 
Jet : 


and for the integers n and 2 restricted to the range Lens A and 


1$ KEE? 
m=] 
cos 235 cos njS = O (1 n) 
d , 
"SD" aa) (2 =n) 
and 
m-.L 
> sin dÉ - Ch Sin dÉ - js = 0 (2 # n) 
j-0 


t 
hol 
d" 
+ 
o 
= 
MEL 
PM 
e 
H 
E 
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The definite integral, which is derived in appendix C, is 


l 
O 
~~“ 
eJ 
"ll. 
B 
Sen 


i | 
J g5( 1,6) Ro (n, d)af 


i 
P 
H 
+ 
o 

E 
IB 
m~ 

e 
u 
B 
c 


where n and 2 are restricted to 2 <n S 
simplification the following equations result: 


MIE 
3 
no 
IA 
e 
IA 

ME 
: 
d 
ky 


For i= 1, 


fn(3) + (27, - 1)£4(2) + 2(Bn - yn)fn(1) = aal ls n5 


3 2 
(46a) 

for 1-22, 
fn) + Qypfn(3) + 2B84f4(2) + (27, - 1)fa(1) = 0 (46b) 


and, for i 2 3, 
fn(i-2) + 2y7pfn(i+1) + 2Bnfnli) + 2ynfn(i-1) + f,(1-2) = 0 (46c) 


Solution of finite-difference equation.- Equation (46c) is the same 
finite-difference equation for which the solution is written in the two 


preceding sections. Substitution of this solution, equation (38), into 
equations (46a) and (46b) gives two simultaneous algebraic equations 
for ajy and Qop; the arbitrary constants. Solution of this system 
yields 


Bn - byn - ll n5 
T BE A Hoc o omm a 
2n z sin 5 " | 


Qin 7 ~ 


(11gl'og =- Pinon an T 55,2) 
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Bn ~ 4yn - 11 nd 


r ee gin See 
"in 5 NE Lar 
en ~ o. T T..0 | . 
inten 7 Tinton a +5 2) 
n5 


The Q's and I's in this case are precisely the same as in the pre- 
ceding case of the distributed perturbation load; Lan is given by equa- 
tion (39a) and Dan by equation (59b). 


With the coefficients fQ(i) known for the shear perturbation load, 
the stringer loads are obtained from equation (40) and the shear flows 
can be found from equations (41) and (hk). For panel row j = O, the shear 
flow equations become 


m or m-1l 
2 2 fa(i-l) - fali) 
T EE (121) 
10 n5 
n- 2L sin = 
2 
and. 
m m-1 i 
2 Or e 
3g, * | fn(1) 29 
doo = — + Ss + ———————— 
n- L sin X SL + A 


When the shear flows in panel row j = O are known, it is simpler to 
compute the -remainder of the shear flows by use of the equations of 
statics rather than equations (41) and (44). In shear panels (0,1) 
and (0,-1) adjacent to the loaded panel, the shear flow is given by 


E 2P]] + QL 
3017 90,-17 9007 — T 
All the other shear flows are found by use of equation (10). If desired, 
the ring bending moments can be found from equations (42) and (45) and 


the ring thrust and transverse shear can be calculated from the formulas 
given in appendix B. 


ha 


gà 


where 
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LIMITING CASE OF RIGID RINCS 


If the ring bending stiffness is allowed to increase indefinitely, 
the rings approach complete rigidity in bending, the parameter C approaches 
zero, and a considerable simplification results. For this limiting case, 
equations (20) for the concentrated perturbation load reduce to 


f,(itl) - 2 a fylt) + fg(1i-1) = 0 (121) (47) 
n 


where 


An = 5B5^ - 1+ cos nô 
B, = 385? + 2(1 - cos nō) 


This can be shown easily by multiplying equations (20) through by C 
and allowing C to approach O. Equation (47) is a second-order finite- 
difference equation with constant coefficients. The same equation, 
together with its general solution, is given in reference 6, page 31. 
The solution compatible with the boundary conditions at infinity can be 
written as 


fn(i) = lte ef (48) 


cosh Ay, SÉ 


and where the upper sign is taken when A, > O and the lower sign when 
Ay, <0. 


Ihe arbitrary constant ayn is determined by evaluating the solution, 
equation (48), for i - O and introducing the value of f,(0) given in 
equation (9). The result is identical to a), given in equation (22) 


z F > 
n= 2 
On (n2 ) 
72 
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Equations (8) and (12), the expressions for stringer loads and shear 
flows, respectively, used before in the case of the concentrated pertur- 
bation load are still valid. The substitution into these expressions of 
the solution (48) with the constant Gr as found above yields the stringer 
loads and shear flows due to a a Lae Peewee load when the 
rings can be considered rigid. = 


For the case of the distributed perturbation load, equations (37) 
reduce in the limit to 


(-An)fn(2) + (2Bn + An)fn(1) = 3B8°Sd, 
fr(itl) - 2 = fr{i) + fQ(1-1) = 0 (1 2 2) 
The arbitrary constant an in the solution (48) is 


GES S 
= m/i + 5, D 
Al +e An + 1) vi A 
For the shear perturbation load, equations (46) reduce to 


(-An)fn(2) + (Bn + An)fn(1) = -6LQigBb^sin © 


fPr(itl) - 2 = fali) +-fnli-l1) = 0 (1 2 2) 


The solution is again equation (48) and a, becomes 


IW 


As 
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CONCLUDING REMARKS 


Solutions have been obtained for the stress distributions caused 
by three types of loading on infinitely long, eircular, semimonocoque 
cylinders with flexible rings. The solutions are adaptable to the con- 
struction of tables of influence coefficients which give stringer loads 
and shear flows in the neighborhood of the perturbation load due to a 
unit magnitude of that load. The perturbation loads are the tools of a 
method described in NACA Technical Note 5200 for the stress analysis of 
circular semimonocoque cylinders with cutouts. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., April 22, 1954. 
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APPENDIX A 
SUMMARY OF SIGNIFICANT EQUATIONS 


The formulas and parameters required for computing the stress dis- 


tribution due to concentrated, distributed, and shear perturbation loads 
are collected in this appendix for reference. 


Stringer Loads 


Concentrated perturbation load (see fig. 1(a)): 


m m-l 
299 


pij = Ap * Boos b > fg(i)cos njb (1 2 0) 


Where P is the applied load. 


— 


Distributed perturbation load (see fig. 1(b)): 


m m-1 
= or —— : 
Pij = d cos jo + > Zoll loop nj8 (1 2 1) 
em m E 
* 
where S9 is the total applied load. S , , 


Shear perturbation load (see fig. 1(c)): 
2 x S 
Pj= 2 f,(i)sin d e is (1 2 1) 


Shear Flows 
Concentrated perturbation load (see fig. 1(a)): 


For the shear panels in panel row j = O, 


— P109 7 P4413,0 
qo e —————— 


2L 


NACA TN 3199 ` 


and, for the remainder of the shear panels, — 


P13 T Pit, 
dij = Z + 04,4-1 (3 


IV 


Distributed perturbation load (see fig. 1(b)): 


For the shear panel (0,0), 


_8 = 2pio 
2r 


doo 


for the remainder of the shear panels in panel row A = O, 


_ Pio 7 Pi+l,O 
uA e 


"nV 
ES 


and, for all other shear panels, 


Pij 7" P111,j 
dij = t+ Oy j-1 (j Z 1) 


Shear perturbation load (see fig. l(c)): 


For the panel &bout which ihe load is applied, 


m m=} 
nm fu 1) " 2Q 


n-2 L sin = En +O, n) 
A 


VN 


doo = AR 


for the remainder of the shear panels in row j = O, 


m or Mel 
2 2 f,(itl) - fy(1) 
dio = EFE (i 
n= 2L sin r3 


V 


1): 


1) 


22 


36 


for the shear panel (0,1), 


301" 9007 e 


for the remainder of the shear panels in panel row j = 


a Pil 7 Pirl,1 
41 ~~; t dio 


and, for all other shear panels, 


Pij ~- E * 04,31 


Evaluation of the Trigonometrie Coefficients fli) 


for Flexible Rings 


Basic parameters: 


Ud 
I 
GJ 
| 


Q 
H 


Auxiliary parameters: 


Bn =3 + 
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P 


(12 1) 


(34 2 2) 


CH 


NACA 
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Discriminating parameter: 


_ 2(Bn - 1) 


Yn 


Trigonometric coefficients: 


fali) = bn" biss?) = Segen (1) | 


Xn 


Yn 


"^n  -*n 


Tal 


bn =- 


Ain(i) = cos ix, 


= 1 


cosh ix, 


Aon(i) = sin ix, 


= sinh iX, 


(D, > 1) 
(D. = 1) 
(Da < 1) 


2T 


(D, > 1) - 


(Dy = 1) 
(D, < 1) 


(Dn > 1) 


(Da < 1) 
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Arbitrary constants for concentrated perturbation load: 


€5 


P 
a + 55,2) 
7 94n + 2(yn +L) P 
Gon m/l-4 8, m) 
C tu) 
where P is the applied load and. 
San = t Asa (3) + 2y ta Aen (2) + (28, e 1) E Asn (2) (s = 1,2) 


Arbitrary constants for distributed perturbation load: 


Bn - Jon - 2 
fon cope Fan 28 
Din ER 
Lil an 7 Tin®en Wi + bn 8) 
By ~ äs - 2 = 
oan = - ————————— ———— 
Hiel an 7 linfen dp + 55,2) 
where S is the total applied load and 
asn = Ee Aen 2) + (27n - 1)tn Asn(2) + 2(Bn - Yn) tn^sn(1) (s = 1,2) 


Den = Kass? + 2y, tn Men (3) + 2Bntn Aen (2) + (27n - 1) pAgn(1) 


(s = 1,2) 


t? 


NACA 


where 
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Arbitrary constants for shear perturbation load: 


Bn - n - 11 
Da — Sin Bo 
3 2 hot, 
Qin = ep. ee 
Wio zo 7 Tinton eh +6, m 
"2a 
B, - ky, = 11 
Se n | n sin È 
E 3 2 Jor, 
Ai oan ~ Tinton Séi + 85,2) 


Q is the applied load per unit length. 


Evaluation of the Trigonometrie Coefficients f,(1) 


for Rigid Rings 
Basic parameter: 


AR 


Baz 
G bt I2 


Auxiliary parameters: 


An = 3B57 - l + cos nô 


B, = 3B5? + 2(1 - cos nd) 


Trigonometric coefficients: 


1 
socie 


29 


Lo 


Arbitrary constant for concentrated perturbation load: 
P 
mfl +65 
C t») 


Arbitrary constant for distributed perturbation load: 


Qn = 


=A m/l +8 m 
S e “+ 1) ( EH 


Ki 


Arbitrary constant for shear perturbation load: 
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APPENDIX B 
BENDING MOMENT, AXIAL THRUST, AND TRANSVERSE SHEAR IN RINGS 


Expressions will be developed for the bending moment, axial thrust, 
and transverse shear in a circular ring under tangential loads such as 
those which arise from the differences in shear flow across a ring in a 
circular semimonocoque cylinder. 


Iwo cases must be considered: One case occurs with the concentrated 
and distributed perturbation loads, where the ring loading is anti- 
symmetric about stringer j = 0. The other case occurs with the shear 
perturbation load, where the ring loading is symmetric about panel 
row j = 0. 


Concentrated and distributed perturbation loads.- For the concentrated 
and distributed perturbation loads, the tangential loading on ring i has 


been written in the form of a finite trigonometric series (see eg. (15)) 


M or ml 
Fes = - "s du sin n(j + Ze (21) 
LIED €r121.] Gs ni 2 


where 


Ayifn(i) 


no 
eL sin 5 


Phi = (n 2 2) 


This ring load has a stepwise variation around the ring, being constant 
between stringers and having jump discontinuities at the stringers. The 
limitation that n 22 ensures that the ring is in equilibrium. 


The procedure will be to expand each term of the series (Bl) in an 
infinite Fourier series in the variable Ø. For each harmonic of the 
Fourier series, that is, for a continuous sinusoidal tangential force 
distribution on the ring, the moment, thrust, and shear in the ring are 
easily found (see ref. 5, p. 55, for example). On the basis of inexten- 
sional deformation and the neglect of transverse shear distortions, the 
results are as follows: If the tangential load on ring i is given by 


Sun cos nf + by, sin nj ` (n 2 2) 
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then the moment, thrust; and shear in this ring sre, respectively, 


Doll «dia ee PE NEN d 
li = en E ————- Sin Y) + b — COS N 
Matt» ni [G2 - 1) ni (m? - 1) 
— R — R 
Tnli,d) = -a&i ——— n sin ng Lk Dat ———— n cos ng (B2) 
n? = i n2 = 1 
Vi (1,8) =a, -—— nø + b. Li ui ng 
n İ p2- 1 ni ne L . 


Figure 5 shows the sign convention used in writing equations (B2). 


Consider, now, one term of the series (Bl). To expand this term 
in & Fourier series, write e 


Dan sin zl + Ze = ER (cr); sin rø | (B3) 


r=n 


where the Los lui S are the Fourier coefficients. It is obvious that- 


the first harmonic which will occur in the Fourier series in equation (B3) 
must be that for which r =n. The other harmonics, then, will be added 
to this to build up the step.shape of the loading function. The con- 
vention for measuring angle Ø in this case is illustrated in figure (a). 
The index j can be thought of as a function of Ø, that is: when 
O<$¢@<8, j=0; when 5< G< 25, j= l; and so forth. 


In order to carry out the expansion, equation (B5) is multiplied 
through by sin ij and integrated from O to 2x 


J=0 "jo 


Ep bni sin aC + Bp sin 29 dg = J” 2 (cr) sin rø sin 1¢ ag 


After integration, the right-hand side of this equation becomes 


(cz) 4% 


D 


e? 
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by virtue of the orthogonality of the trigonometric functions. The left- 
hand side becomes 


2 sin = a 1 1 
EAE RE bui sin dÉ is z) sin dÉ + Lp 


Es 


on carrying out the integration. From reference 4 it can be shown that 


m=] = | tin 
1 L -m m 
sin n(j + =)S sin 2[j + =\8 = 2 {(-1) Jj.4 - (-1) J 
a Oh ieee ee | i stad "e 


where Jy, =1 if r is an integer, and J, =0 if r is not an integer. 


Thus the Fourier coefficients are given by 


io 15 i-n ttn 
m es m m 
(ez) = F Pmt —4-—j(-3) pe C1) Tum 
m m 


The nth term of the tangential loading on the ring is 


Se i-n itn i 15 
m 
bni sin n() + Ab = B bni 2 -D E Za CD)" ona sin 1 
=n m m 


By use of the properties of J, this summation can be rewritten 


sin(rm + n) 


li. -m S T 
bni Sin dÉ + Sie = = bni Es (-1) OU ES sin(rm + n)¢@ - 


co in - yo 
S- Qr LUE sin(rm - n)f (85) 


r=l 
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On expansion by the sum and difference formulas of trigonometry and with 
the use of the fact that—m5 = 2x, it is found that- 


(~1)"sin m 


sin(rm + n)? 


(B6) 


5 r+l no 
sin(rm - n)2 (-1) '^sin = 


When equations (B6) are substituted into equation (B5), the following 
relationship results: ` 


SE ht . 2j FEMME" 2 sin(rm + n)ó E = sin(rm - n}¢ 
ni 2 TÉ ni 2 Es m+n v rm-n 
CO LI 
m nS sin(rm + n)¢ 
xe 2 a. m+n (BT 


From the first of equations (B2) it is seen that if-the tangential 
loading on the ring is given by the right-hand side of equation (BT) 
then the bending moment in that ring is 


M,(i,g) = R? S ba; sin 2 Hi (n,@) (B8) 
where 


oO 


WE EE en 
a" (rm + SE en + n)* - i| 


Equetion (B8) gives the bending moment in & ring which carries & 
tangential load distributed according to one-term of the series of 
equation (B1). When the ring is loaded by the sum of such stepwise 
terms, as in equation (Bl), then the moment is given by a sum of terms 
like (B8). So the bending moment in ring i is 


"së Eon sin È wy (n,9) (29) 


vd 


c? 


') 
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For completeness, the expressions for axial thrust and transverse shear 
can be written in a similar manner - 


m m= 1 
peo 
T(1,9)- > R B bai sin È K1(n,) 
n=2 
B or Bl 
V(i,g) = R Bt sin È 1y(n,9) 


where 


TP —— cos(rm + n) 
Ping Ke -1 


CO 


= sin(rm + n)9 


(n,g) = ————————— 
D | r=- (rm + 2) [(rm + n)? - 1] 


Shear perturbation load.- In the case of the shear perturbation 
load, the tangential loading on ring i is given by the finite trigo- 
nometric series 


m m-1 
ie 


2 
Fij = j- -g= 2, ent cos nj — (B10) 


where 


_ Aiifnli) 


no 
eL sin 5 


(n 2 2) 


Equation (B10) can be treated in a manner analogous to the handling of 
equation (Bl). That is, each term of the series in equation (B10) can 
be expanded in a Fourier series. Then the moment, thrust, and shear in 
the ring are written immediately. 
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Analogous to equation (B5), write-— 


La 


co 


Soit cos një = 2. Los Jeu oe rø EN (B11) 


r-u 


where, now, the angle d is as shown in figure h(b). If both sides of equa- 
tion (Bll) are multiplied by cos 10 and integrated from O to 2x, there 
results for the Fourier coefficients: 


2 sin © m-l 
(ez)ni = — ani > cos nj8 cos 1j5 
j=0 ' 


It can be shown (see ref. 4) that 


cos nj5 cos ijo = 2 daat d 
m m 


so the nth term of the tangential loading on the ring is 


Go sin 5 
ani cos njẹ = = ani > din + dän Pe cos 2¢ (812) 
LO Um m | | " 


This summation becomes 


Qo 


Se E SS > tt cos(rm + n)¢ 


ani cos nj5 = 
rz.co rm + n 


AJB 


which corresponds to equation (B7). Then the benđing moment is 


m or t-l . 
2 2m nB & f. 
M(1,9) = - > R^ 7 ani sin 5 Hp(n,@) (B13) 
ne Te. 
Similarly, thrust and shear are 
mor Bat 


2 | | . 
T(1,9) = - 2. R = api sin 2 Ko(n,9) 
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" 
] and 
5 or Bet 
v(i, Ø) = "3 R x ani sin 2 Lo(n,9) 
where 

Staf >. (1 si 
2. (rm + 2)? |(rm + n)? - 1] 

blue a a 

y=- (rm + wo - 1 

TA F S cae cos(rm + n)Ø 


(rm + n) |(rm + n)? - 1| 
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APPENDIX C 
EVALUATION OF DEFINITE INTEGRALS 


In order to minimize the stress energy it is necessary to investigate 
the following definite integrals: 


L 81(,9)8 (2,9)ag = I 2. Den cos(rm + n)g ES EE up qp 


Ee EEN 


(c1) 
and 


2 ent e co 
I i Ho(n,$)Ho(2,9)ag ei X (-1) Den sin(rm + n)¢ S (-1)*°D,; sin(sm + 119 ad 
D'So SES ent 


(c2) 


where 


l 
een e oe 1] 


and where integers n and 2 are limited to the following ranges: 


Drn = 


<£ 


HA 
H 


e 


e 


MIB mg 


n^ 
e 
WA 


Consider the relation (C1). The right-hand side can be written 


ex DQ 
z P3 2. Halle 1 cos(rm +n - sm- 1) ag +f cos(rm +n + 
H O 


ne) | ` 


y--90 g= 


co GO 
sm + 1)9 d = > > Deng 1 (Bremen, amt? = Srmin, -sm- 2) 2x 
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Now, by virtue of the limited range of the integers n and 1, the following 


relations can be written: 


e? 


Doan, emt) x E. T ten ^ Osr? in 
3 mb ` "-——— 
n 
Srmin,-sm-2 = GC ln = Ss ,-r-157,@ Ba 


Thus, when 2$ De 2: equation (Cl) yields 


l 
O 


ex 
J, Hy (n,Ø)E1(2,ø)ag = 


S 2 
Drn i = Snr 
Y= 


If n= > the following equation is obtained: 


- I EAG Aag = > gs "Top EE 


Since 
E 
(m-at Boma a 
E es 
SCHEER 
^ = D | 


(2 # n) 
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it is Pound that when n= 


MIR 


Ge e 
[ Sg 8, 9 )ag =2 > d vn = Së 
O Y—-co 22 2 


To summarize, then, 


ll 
© 
~ 
e 
Be 
H 
ath 


2x 
D H1(n,2)11(1,9)ag 


bat? $ 5,1) (2 =n) 


Consider the relation (C2). It is handled in a manner analogous to 
the treatment of (Cl). Equation (C2) can be written 


NE T 2. > (-1)"Dry(-1) "Dg 13(Srmin, sm? D Boun, -em- 1 27 


oO 


S > KR LATE i (55:81 - RTE 


Y==400 £g--oo 


For DS n«z 


li 
O 
7 
~ 
wW 
B 
— 


en 
J, Ho(n,9)B5(1,9)ad 


Il 
Co 
^ 
A 
il 
d 
~ 
e 
il 
CO 
Kg 


LA 


ei 


Le 
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[, ota - 


M 
[S1 


2 r -f= 
dÉ zo (a) Dyn D d 


-r-1, 


I 
Me 
| oz | 
NIB 
nó 
| 
T 
= 
| 
d 
H 
à AA 
a 


l 

no 

Jg 
B 
a 


Il 
tà 
a 


Thus the definite integral (C2) gives precisely the same result as (C1) 


ex 
JL Taepa = o dP 
= "et + oe (2 = n) 
The sum 
Sn = Din = 
RN Sp (rm + St? D + i) S i 


can be expressed in closed form with the aid of formula 6.495 ; number 2, 
reference 7. The result is 


_ Si 2 + cos nô 5? B^ cos nb cos B - 1 | 5 5 sin 65 


+ Ee 
12 (1 =- cos nb)? 1 - cos nô 4 (cos nB - cos 8)* 4 cos n - cos 8 


However, the series form of Sp, because of its rapid convergence, may 
be more convenient than the closed form for use in computation. 
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Ring Ring 


(c) Shear. 


Figure 1 .- Perturbation loads. 


Ring i-2 i- | i i+] i+2 


Figure 2.- Portion of typical cylinder. 
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Detail of ring element. 


(a) For concentrated and (b) For shear perturbation load. 


distributed perturbation loads. 


4% 


Figure 4.- Conventions for angular coordinate œ. 
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